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Abstract. Two-level electron-phonon systems with reflection symmetry linearly coupled to one or two 
phonon modes (exciton and Ecg>(6i + 62) Jahn- Teller model) exhibit strong enhancement of quantum 
fluctuations of the phonon coordinates and momenta due to the complex interplay of quantum fluctuations 
and nonlinearities inherent to the models. We show that for the complex correlated quantum fluctuations 
of the anisotropic two-level systems the Shannon entropies of phonon coordinate and momentum and their 
sum yield their proper global description. On the other hand, the variance measures of the Heisenberg 
uncertainties suffer from several shortcomings to provide proper description of the fluctuations. Wave 
functions, related entropies and variances were determined by direct numerical simulations. Illustrative 
variational calculations were performed to demonstrate the effect on an analytically tractable exciton 
model. 

PACS. 71.38.-k Polarons and electron-phonon interactions - 63.70.+h Statistical mechanics of lattice 
vibrations 

1 Introduction proaches: Feinberg et al ^I] m a model of an exciton, 

Borghi et al m the Holstein-Hubbard two-level model, 
In spite of long-term research, various aspects of the the- Morawitz et al [Hj in the Holstein-Peierls model, 
ory of polarons with Holstein coupling in two-level local Recently, we have investigated numerically and vari- 
and lattice models related with quantum fluctuations and ationally the interplay of quantum fluctuations and non- 
phase transitions belong to systematically studied topics linearity inherent to the two-level models in the ground 
in the current literature up to the present time 0-0- state of the E®(i>x + b 2 ) molecular (local) JT model |H1 
The class of two- level electron-phonon models with lin- and the respective lattice JT model • Corresponding 
ear coupling to one or two phonon modes represent several Hamiltonians contain a hidden nonlinearity due to the re- 
important physical systems: excitons, dimers, Jahn- Teller flection symmetry (the nonlinearity appears explicitly un- 
(JT) systems (rotationally symmetric E®e and reflection der an appropriate unitary transformation; See Section 
symmetric E®(&i + b 2 ) JT model with broken rotational 2 and discussion in our recent papers 53,^3). There is 
symmetry). Especially, interest in the literature is grow- to elucidate the terminology related to quantum fluctua- 
ing in the JT models [21,0; There is the experimental tions and the nonlinearity parameters used there and 
evidence of related structural phase transition in some in the present paper: quantum fluctuations are measured 
spatially anisotropic complex structures (perovskites, ful- by the ratio of the phonon frequency and the classical pa- 
lerides, manganites 0,13,111: EH)- rameter of electron-phonon interaction, Q/ot = l/y/2jl, 

In two-level electron-phonon models strong enhance- while the nonlinearity parameter is the quantum tunncl- 
ment of fluctuations of phonon coordinates and momenta ing strength between the levels, /a = x- I n the plane 
as well as of their product in certain ranges of model pa- /i, x (p — a 2 /2Q 2 ) there occur quantum fluctuations ~ Q 
rameters were reported by several authors from both nu- and <~ x- The ground state in the phase plane fi, x ex- 
merical simulations and analytical (e.g. variational) ap- hibits regions of dominance of either selftrapping (clas- 
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sical) (/i > 1, x < 1) or tunneling (quantum) parame- 
ters ((/, < 1, x > 1) and regions where the classical and 
quantum regions mix together (fj, < 1, % < 1 or fi > 1, 
X > 1). For fi < 1, formation of both the selftrapping and 
tunneling "phases" is suppressed by the quantum fluctu- 
ations oc J?. The parameter of nonlinearity x and phonon 
frequency Q are both of quantum origin so that global 
quantum fluctuations mix the fluctuations due to nonlin- 
earity close to the border of two regions at % ~ 1 and the 
quantum fluctuations ~ S7. 

Numerical approach to the two-level exciton and JT 
models in consideration in comparison with variational 
approaches bring understanding of the nature of mixing 
of quantum fluctuations with the nonlinearities: since the 
coherent phonon subsystem does not conserve the number 
of phonons the upper level participates in the distribution 
of phonons even in the ground state. This manifests itself 
by the appearance of additional reflective extrema of the 
ground state wave function |16j.|17j. The distance between 
its extrema is related to the displacement of the coherent 
phonons and is involved into the dispersion measure of the 
wave function and of its Fourier transform |18j . |19j . |20|. 
This displacement, i.e. the measure of the polaron selflo- 
calization, though suppressed by quantum fluctuations, is 
of classical origin. This is the source of difficulties with the 
justification of the moment (variance) characterization of 
fluctuations in the case of the wave function with the pres- 
ence of the additional reflection maximum. As it is known, 
in the case of multipeak wave functions, the variances for 
the coordinate and/or momentum of related wave func- 
tions, as well as for other non-commuting observables do 
not stand for appropriate uncertainty measures [TS| ~ |22J > 
since these uncertainty measures strongly involve distance 
between the peaks. 

Another serious shortcoming arises when we switch 
from one- to two- or more phonon system, like Jahn- Teller 
systems are (generally - to a system with several non- 
independent variables). In this case the moment-related 
uncertainty measures are not good even for one-peak dis- 
tribution, since the width of the distribution given by the 
variances is a width in a particular direction chosen arbi- 
trarily, and cannot stand for actual "width" of a distribu- 
tion. 

In these cases the Shannon entropies of the probabil- 
ity density functions assigned to such wave functions were 
used as alternative uncertainty measures of the conjugated 
coordinates. The "entropies" of any probabilistic distribu- 
tion are well-known to eliminate effectively both the dis- 
tance between peaks in the multipeak distributions and 
the said "anisotropy" of the distribution in the space of 
random variables since they contain in their expressions 
only functions of probability distributions, rather than val- 
ues of a random trial, as in momentum measures. 

For references, we shall briefly summarize necessary 
notions: in the probability theory there are two main types 
of the integral uncertainty measure assigned to an observ- 
able, the moment and entropic ones J25J- Due to existence 
of these uncertainty measures two types of uncertainty 
relations (UR) for two non-commuting observables can 



be introduced: (i) the Heisenberg (variance) UR and (ii) 
the entropic UR. While the moment (variance) UR is ex- 
pressed as the product of the variances of two noncommut- 
ing observables A and B, the entropic UR is given by the 
sum of their Shannon entropies. For a continuous observ- 
able A with the density of probability distribution p{x) 
the (differential) entropy of A is defined as [21] 

H c = - J p(x) \ogp(x)dx . 

The entropic uncertainty relation for the coordinate 
and momentum of a quantum system described by its nor- 
malized function ip{x) is represented by the inequality |25j . 

S x + S p > S xp , (I) 

where S x and S p are the entropies of its coordinate and 
momentum probability distributions 

oo 

S x = - J |V(z)| 2 log|V(z)| 2 cfe (2) 

■ — OO 

and 

oo 

S P = - I \^p)\ 2 \og\^(p)\ 2 dp, (3) 

— oo 

respectively, and <p(p) is the Fourier transform of the wave 
function ip{ x ) and S xp represents the lower bound of the 
right-hand side of the inequality £[J . 

It has been found by Bialynicki-Birula and Mycielski 
[2E| that the lower bound S xp = 1 + log7r for the har- 
monic oscillator represents the minimal lower bound for 
any quantum system with one pair of non-commuting ob- 
servables. Therefore, the sum of coordinate and momen- 
tum entropies of the arbitrary quantum system is (h = 1) 

S x + S p > 1 + logTT. (4) 

This relation is an entropic counterpart of the famous 
Heisenberg principle formulated by means of variances: 
AxAp > 1/2. 

We note that the entropy as a measure of uncertainty 
and the entropic uncertainty relations are widely used 
in quantum optics [22, ESI being useful for systems with 
more complex structure of photon coherent states (e.g., 
Schrodinger cat coherent states). 

In Section 2, two standard two-level electron-phonon 
models with linear electron-phonon coupling are described 
as prototype models for which the entropic uncertainty 
principle represents the adequate measure of the fluctua- 
tions. The nonlinearity hidden in the (initially linear) two- 
level phonon Hamiltonians is revealed after appropriate 
diagonalization of the problem with respect to the elec- 
tronic subspace. 

In Section 3 we determine by direct numerical simula- 
tions the wave functions of one-phonon exciton model for 
a range of model parameters. From those, we determine 
and compare the quantum fluctuations as functions of the 
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effective interaction fi described by means of the Shannon 
entropies of the phonon coordinate and momentum on one 
hand and by means of their variances on the other hand. 
Comparison of the sum of the coordinate and momentum 
Shannon entropies and the corresponding product of the 
variances clearly puts forward the former as more ade- 
quate measure of the quantum fluctuations in certain pa- 
rameter region than the latter one. 

Our approach is consequently based on numerical sim- 
ulations, because variational treatments based on squeezed 
coherent states or their linear combinations were generally 
found to underestimate quantum fluctuations : (i) they un- 
derestimate quantum fluctuations originating from finite 
phonon frequency |2Z],|2E1 f° r weak couplings; (ii) The 
fluctuations due to the nonlinearity in the crossover re- 
gion are strongly coupled with the fluctuations (i). There- 
fore, in the region close to the crossover there appear dis- 
continuities which are artefacts of variational approaches 
(similar to those of the adiabatic approximation). 

In Subsection 3.1 leaning upon variational treatment of 
the exciton model we present results of simple analytical 
estimations of quantum fluctuations (both entropic un- 
certainties and variances) to illustrate the failure of vari- 
ational methods. However, in a small range of the phase 
plane we can identify the effect of the classical displace- 
ment in the Heisenberg variances and their product by 
means of variational method. 

Similar investigations for the two-mode phonon E (8 
(b\ + 62) Jahn- Teller model (which differs from the model 
of Section 3 by the presence of phonon assistance in the 
tunneling term and, consequently, the existence of addi- 
tional mode correlations effects) are presented in Section 
4. 



2 Model Hamiltonians 

For the model of interest let us consider the two-level 
reflection symmetric electron-phonon system with linear 
(Holstein) coupling to one or two phonon modes 

JV=1,2 

H= n(b\bi + 1/2) + a{b\+b x )a z - [3Aa x , (5) 

where the Pauli matrices a x , a z represent electron den- 
sity operators in the pseudospin notation, <j z = -(0^02 — 

c\ci), a x =J>(c\c2 + c\ci) with / = c\c\ + cj>C2 as the unit 
operator J^j. 

For the case of coupling to one phonon mode, N = 1, 
we take A = 1 and two-level Hamiltonian (J5J represents 
an exciton or a dimer with tunneling between the levels 
with the tunneling amplitude /3. The prototype model has 
been handled both variationally and numerically by many 
authors in its local and extended version ^3>GB1jGZ1> P^l - 
In the case of the coupling with two phonon modes, N = 
2, the second term with A = b\ + & 2 represents phonon- 
assisted tunneling between the levels (the flip-flop rate 
is proportional to the value of phonon-2 coordinate, i.e. 



an additional "phonon pumping" occurs in distinction to 
the one-phonon case). The special case a = (3 represents 
standard E ® e Jahn- Teller Hamiltonian with rotational 
symmetry [HJ , [HI , , EH^ESj- The phonon 1-mode is 
antisymmetric and the 2-mode is symmetric against the 
reflection. The case a ^ f3 is the reflection symmetric 
E®(bi+b2) Jahn- Teller Hamiltonian recently investigated 
variationally and numerically |14 ) .jl5 ) . 

The reflection symmetry of the Hamiltonians j£| in- 
volves a nonlinearity which reveals itself explicitly after 
performing Fulton-Gouterman unitary transformation (Ful- 
ton et al glj, Shore et al QU): 

N=l,2 

H FG = UHU- 1 = n ( b th + 1/2) + a(b\ + h)I 
i=i 

TPAexp{iirb{bx),(6) 

where U = ( i'^qj j G = exp{iirb\bi) is the re- 
flection operator in the phonon space, G(b\ +61) = — (b\ + 
bi)G, which performs virtual coupling of the levels by 
phonons 1 mediating the electron (Rabi) oscillations be- 
tween them. 

The transformation (jfJJl diagonalizes (J5J in the elec- 
tronic subspace yielding however a strong nonlinearity in 
the phonon subspace (term containing (3) which otherwise 
was hidden in the initial Hamiltonian JSJ). The unitary 
transformation left us with purely phonon Hamiltonian 
© while electrons were excluded. All further consider- 
ation will be performed for the ground state which for 
both models is given by the upper sign in the transformed 
Hamiltonian |JBJ. 

The competition of quantum fluctuations (J?), the self- 
localization (a) and the phonon assisted tunneling (/3) in 
© determines regions of dominance of said effects for dif- 
ferent sets of the parameters. In the next sections, we shall 
analyze these regions by numerical simulations of the mod- 
els ©. 



3 Exciton: coupling to a single phonon mode 

The main purpose of this and of the next section is to show 
that for the considered prototype models the entropic un- 
certainties and the entropic uncertainty principle repre- 
sent adequate measures of quantum fluctuations instead 
of the moment (variance) Heisenberg ones. 

The exciton model, including respective wave func- 
tions, has been previously intensively investigated by many 
authors; Here we shall confine ourselves to essential points 
needed for our purposes. 

The results of the numerical simulations of the ground 
state wave functions of the model © for A = 1, i.e. of 
one-mode diagonal coupling case, are shown in Fig. 1. The 
parameter of the effective interaction fj, = a 2 /2Q 2 is the 
measure of the competition between the classical polaron 
selflocalization due to the interaction of the energy a 2 jlfl 



4 



E.Majermkova et al.: Entropic uncertainty measure 



(e.g., Holstein [jS]) and of the quantum fluctuations of the 
energy fl. Parameter (3/ f2 is the measure of the compe- 
tition between the quantum term of the tunneling [3 and 
quantum fluctuations Q. Wave functions shown in Figs. 
1 illustrate the effects of the interplay between these pa- 
rameters affecting the population of the phonons on both 
levels. One can see that the most prominent deviations 
from the one-peak wave function occur for large For 
moderate /i and (3, (Q — 1) the wavefunctions also tend 
to exhibit variations from almost- Gaussian shapes. These 
features are consistent with non-conservation of the num- 
ber of the coherent phonons involved. They inspired the 
variational approach based on linear combination via a 
variational parameter of two harmonic oscillators related 
to both levels ^H]>GZ]- Numerical evaluation of the ground 
state (Fig. 2) illustrates smooth decrease of the energy as 
function of fi indicating a smeared crossover region. The 
sharp transition line between two regions resulting from 
the competition of two interactions (fi and (3) is known to 
be an artefact of the adiabatic approximation and often of 
inadequate variational approaches. In the quantum treat- 
ment the crossover line is smoothed over the region of a 
width w Q of the phonon energy. 

These regions are analogous to the "selfrapping dom- 
inated" and the "tunneling dominated" regions for the 
E®{b\ + 6 2 ) model UH,[IS| (see the next Section). In the 
weak coupling region, fi < 1, the formation of the "ordered 
phases" (either selftrapping- or tunneling- dominated) is 
suppressed by the quantum fluctuations ~ Q as well as in 
the model of Section 4. 

Corresponding to the above described reasons for qual- 
ifying the quantum deviation from the harmonic oscillator 
(Gaussian) behaviour we evaluated numerically (Fig. 3a) 
the left-hand side of the Shannon uncertainty relation, the 
sum of entropies Sq + Sp (@J (here, Q and P are the co- 
ordinate and momentum of the phonon). The product of 
variances y/ AQ 2 AP 2 containing explicitly the displace- 
ment (distance between the peaks) is shown on Fig. 3b. 
Both approaches, the variance and the Shannon entropy 
ones, can be compared to illustrate their different physical 
content. Close to the line of equal effective coupling and 
the tunneling strength 2/i = /3 the extremum of the sum of 
Shannon entropies manifests itself markedly. The pictures 
of related wave functions attribute this extremum by the 
maximal overlapping (tunneling) of the contributions of 
both levels (Fig. lb). With growing polaron selflocaliza- 
tion \x the quantum fluctuations are suppressed and both 
oscillator parts tend to separate, see Fig. 3a for large fx. 
The harmonic oscillator value 1 + log7r — 2.14473 is re- 
covered in both regions around the crossover region: in 
the limits of large /i or, on the contrary, for large /3 and 
small or moderate fi (less than some critical value on the 
crossover line). 

The product of variances AQ 2 AP 2 at Fig. 3b shows 
similar behaviour except for large /i, where it therefore 
accounts for the growing distance of the peaks of the wave 
functions with /i, Fig. Id. This distance as a function of a 
classical parameter n is a measure of the nonlinearity since 
it maps the position of the reflection peak. In the limit of 



small /i, y AQ 2 AP 2 tends to the harmonic oscillator value 
0.5. 

The range of maximum fluctuations in the plane /i, f3 
(Fig. 3) separates the quantum fluctuation dominated re- 
gion of the harmonic oscillator Q (2/i < (3) and the self- 
trapping dominated region (2fi > (3). Although both un- 
certainty measures yield similar behaviour in the crossover 
range showing the common increase of fluctuations, and 
in the tunnelling region where the wavefunctions indeed 
converge to an oscillator, the quantitative information be- 
yond the crossover region is different in the selftrapping 
region. Enhancement of fluctuations reported by the vari- 
ance uncertainty picture in the selftrapping region point 
merely on the classical contribution related to the dis- 
placement so that it fails to play a role of a measure of 
quantum fluctuations. This shortcoming is effectively elim- 
inated by the Shannon entropic relations which thus help 
to extract purely quantum effects manifested in our case 
in the strongly non-Gaussian behaviour of single peaks. 

Figs. 3c, d compare crossections of both Figs. 3a, b 
as functions of /i for different (3. The entropies and vari- 
ances plotted on one graph (both shifted by the oscil- 
latory values of corresponding quantities) illustrate the 
difference of two approaches: the Heisenberg product of 
variances with increasing fi (upper curves in Fig. 3c, d) is 
growing while the sum of entropies tends to the oscilla- 
tor value as expected. (For quantitative comparison we 
should plot rather exponents of the sum of entropies than 
the entropies; this becomes clear if we remind of the fa- 
mous Einstein relation between entropies and fluctuations 
exp(S') ~ AX 2 , but the difference between them is almost 
inappreciable). 

In the next Subsection we shall illustrate these consid- 
erations by means of analytical calculations of both mea- 
sures based upon a variational fitting of wavefunctions. 

The enhancement of quantum fluctuations for the ex- 
citon model have been found by Feinberg, et al. jl lj by 
numerical calculation of the moment (variance) uncertain- 
ties and their product as functions of parameters a 2 = 2/i 
and A = 2/i//3. For given a, their dependence on A means 
in our notations the dependence on 1/(3. Let us note that 
these fluctuations are related to the competition of the 
parameters /i and of the nonlinearity parameter (3, i.e. re- 
lated to a crossover between the selftrapping dominated 
and the tunneling dominated regions. However, the un- 
desirable presence of the classical nonlinearity manifests 
itself in the dependence of \J AQ 2 AP 2 on fi due to the 
competition between the parameters a and Q. Namely, it 
becomes obvious when the product of variances is depicted 
as a function of a 2 = 2fi as in Fig 3. 

3.1 Analytical illustration: variational approach 

Though the variational approaches based on squeezed co- 
herent states are not suitable for description of quantum 
fluctuations we will apply them for evaluations of fluctu- 
ations of a simple exciton model. We can use them for 
demonstration of the effect of interest (i.e. of the involve- 
ment of the displacement in the moment uncertainties 
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product in contrast to its entropic counterpart). We will 
use variational approach with two squeezed coherent har- 
monic oscillators linearly combined by a variational pa- 
rameter (112(1 in accordance with the concept initiated by 
Shore and Sander [TS|. However, this illustrative analyt- 
ics is valid only if the maxima of the wave function are 
well separated from each other (small overlapping), i.e. 
for strong coupling (large fj,). This is just the region of the 
most prominent difference of both said concepts of the 
uncertainty measure. 

Phonon variational wave function for the two-level model 
with one phonon mode (A = 1) as solution to the 
Fulton-Gouterman equation [17 ) . |34 | . |2H | 

H { /1^ = [Q (fit J + l/2)+a(b t +b)-p(3G)]^ p) = 



(7) 
p = ±l 

for the ground state (p = 1) can be chosen in the form 

mm- 

|^M)=-L(l + ^G)^ 1 )), (8) 

where G — exp(i7r6^ b) is reflection operator in the phonon 
space, G|0«) = |^ (2) )G, and 



|^ 2 )( 7 ,r))= D(± 7 )S(r)|0>. 



(9) 



The indices 1(2) pertain to the lower (upper) level. The 
functions D(-j) and S(r) are generators of displacement 
and squeezing depending on the variational parameters: 



D( 7 ) = exp[ 7 (fit - b)}, 
S(r) = exp[r(6 t2 - b 2 )] 



(10) 



acting on the phonon vacuum |0). The normalization con- 
stant in JSJ is C = 1 + rj 2 + 2?yexp[— 2 7 2 exp (— 4r)]. Wave 
function (f^f>— (fTTTf) represents two squeezed displaced har- 
monic oscillators combined by a variational parameter rj 
(compare with Fig.l). 

Variational parameters of the displacement 7 , squeez- 
ing r and the parameter of admixture of the reflection part 
of the upper level rj can be determined by minimalization 
of the Hamiltonian Q averaged over JHJ : 



(if) = - cosh4r + 

-^ 7 2 (1 + rf - 2r, exp (-8r) exp (-2 7 2 )) + 

§(l-^)7-§((l + ^ 2 )exp(-2 7 2 ) + 2r;), 

7 = 7 exp(— 2r). (11) 

(The parameters are scaled so that Q = 1). 

It is worth noting that variational approaches of vary- 
ing degree of reliability combined with unitary transforma- 
tions are widely used for electron-phonon systems |11|— [T7]. 
|22| , [2] , [2U i EHU ■ Comparison of different variational Ansat 
for one-phonon two-level system was performed e.g. by 
Shore et al [TB], Sonnek et al [T7j and for E®(6i + b 2 ) 
model by the present authors |14| . 



Let us evaluate variances of the phonon coordinate 
Q = (fit + b)/y/2 and momentum P = (b f - b)/iy/2, 
(AQ) 2 = (Q 2 ) - (Q) 2 , (AP) 2 = (P 2 ) - (P) 2 , where we 
average over the states (l%))- (|10f) . 

We get 

fAn ,2 1 , A s(, ^ a ~2 2?7+(l + ^)cxp(-2 7 2 ) 
\ A Q) = 77 exp(4r) 1 + 877-7 



(1 + 772 + 277 exp(-2 7 2 )) 2 
(12) 

cxp(-2 7 2 ) 



(AP) 2 = - ex P (-4r) 1 - S^ 2 — — 

2 \ 1 + rf + 27yexp(— 27^) 



by 



(13) 

For large 7, the product of variances can be estimated 



(AQ) 2 (AP) 2 



IQrf 



(14) 



(1 + ?7 2 

where 7 is defined by {ITJ . From {T2J) , iJSJ and JTU it 

is obvious that the anomalous enhancement of fluctuations 
is due to the contribution of the classical displacement 
(reduced by squeezing) 7 and is mediated by the reflection 
parameter 77. If 77 = 0, we are left with a single harmonic 
oscillator as expected. 

It is easy to perform simple illustrative analytical esti- 
mations for the values of variational parameters as func- 
tions of the model parameters. Assuming in (|llfl r, rj small, 
we get an approximate equation for 7 setting d (H) /dj — 
(compare similar procedure in the recent paper |14|). 



7(l + 2/3exp(-27 2 )) = -a, 



(15) 



(from similar considerations the expression for 77 can also 
be found). 

Fig. 4 illustrates results of the variational approach for 
the Heisenberg product of variances: in comparison with 
Fig. 3b, there occurs a complete suppression of fluctua- 
tions in the tunneling region ~ 2fi < [3. Instead of the 
smooth crossover there occurs sharp discontinuity which 
is an artefact because of underestimation of the fluctua- 
tions ~ fl similar to that of adiabatic approximation. In 
the selftrapping region, ~ 2/i > 0, the increase is evi- 
dently caused by the classical contribution because quan- 
tum fluctuations tend to disappear in the classical limit 
of the strong coupling. 

To illustrate the advantages of the entropic uncertainty 
relation over the moment Heisenberg ones for our case we 
calculate approximately the expressions for the Shannon 
entropies of coordinate and momentum for the ansatz JHJ) , 
(|2J). This can be easily done analytically if 77 <C 1 and 
7 3> exp(2r) (last condition meaning that two peaks of 
the wavefunction are well separated and almost do not 
overlap). 

For the entropies Sq and S P ©-© with ip)-lJT71Jl we 

get 

Sq = i(l+log^)+log(l+77 2 )+2r--^log(77 2 )+0( e , 77 3 ), 
2 1 + r/ z 

(16) 
>(-2 
(17) 



zes 2 
S P = l(l+lo g 7r)-2r- ^ +0(e, r, 3 ), e = exp(-27 2 ). 
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From (HI and ((HJl 

Sq + S P = 1 + logvr + log(l + if) - / - 

(1 + 7^)^ 

--^log(ry 2 ) + 0( £l?7 3 ). (18) 
1 + rj z 

For rj — 0, Eq. I|18[) reduces to the single oscillator value 
1 + log 7r as expected. Contrasting l|18[> with (|14fl we see 
that entropy uncertainty relations are weakly dependent of 
the displacement 7 and contain as the main contribution 
owing to the parameter r\ originating from the nonlinear 
effect due to coupling between levels. 

4 Interplay between quantum fluctuations 
and nonlinearity in E <g> (bi + b 2 ) Jahn-Teller 
model. 

Quantum ground state of E <g> (b\ + 62) Jahn-Teller reflec- 
tion symmetric model (a degenerate electron level coupled 
with two phonon modes, one symmetric and one antisym- 
metric against the reflection) has been investigated in our 
recent papers for the one site an d lattice ^Hl case. Let 
us expound the main results relevant to clarify the origin 
of quantum fluctuations of interest. 

Exact numerical simulations of the solution to the Hamil- 
tonian JBJ with A = b\ + b 2 yield the ground state wave 
functions depicted in Figs. 5 a, b, c (in the coordinate 
representation in the space (Qi (8 Q2)), in terms of the 
parameters ^ = a 2 /2Sl 2 and x = P/ a f° r three char- 
acteristic regions |14| . The two-phonon ground state wave 
functions for the case of phonon-assisted tunneling exhibit 
evident mixing of the phonon- 1 wave functions related to 
the lower and excited levels at x < 1) Fig- 5a. Each of 
the levels refers to one of two competing minima of the 
effective potential composed of a- and /3-components, as 
described below, Fig. 6. The mixing due to the two mode 
correlation is pronounced most effectively in the region 
of the dominant quantum fluctuations SI, /J, < l,x — 1 
|14|. At x > 1, the prominent peak refers to a harmonic 
oscillator of the dominant potential well related to the 
selflocalized state of an electron oscillating between the 
levels (Fig. 5c). Because of the phonon-2 assistance there 
occurs a smeared continuous crossover from the regime 
of selflocalization (Fig. 5a) towards the tunnelling regime 
(Fig. 5c) through the intermediate picture close to the 
E<g>e Jahn-Teller case (Fig. 5b). 

In spite of the limitations of variational approaches 
(their failure close to the crossover between two regimes, 
underestimation of quantum fluctuations) they provide 
useful insight into the behaviour of the ground state. The 
complex interplay of the nonlinear and quantum effects 
has been analyzed by numerical simulations and compared 
with results of various variational treatments rjl] . We note 
that the variational wave functions with the admixture of 
the excited symmetric phonon mode (phonons-2) we pro- 
posed recently for the ground state showed significant 
improvement of the agreement with numerical simulation 



results especially for strong e-ph couplings fi. The topol- 
ogy of the effective variational potential, Fig. 6, (Hamilto- 
nian © averaged over trial functions depending on a set of 
variational parameters) is controlled by several model pa- 
rameters; it plausibly can acquire two or more competing 
minima referring to the ground state with a possible ad- 
mixture of the side minimum referring to the excited state 
|14| . As a result, two regions of the ground state appear 
according to which one of two local minima of the po- 
tential dominates; generically two regions are recognized 
- with either dominating selftrapping (x < 1) or tunnel- 
ing (x > !)■ T ne existence of selftrapping dominated vs 
tunneling dominated regions results from complex compe- 
tition of two nonlinearly coupled coherent phonon modes 
(/3-term in ©). Let us note, that the order parameters for 
the selftrapping dominated "phase" is the displacement 
(b\ + 61) ~ Q\ and for the tunneling dominated "phase" 
{b\ + b 2 ) ~ Q2 - These features are analogous for both local 
[j"4"| and a generalized lattice ^H] E®(6j + 62) JT model. 

The additional phonon-2 assistance of the tunneling 
with non-conservation of the number of the phonons 2 
implies their selfconsistent behaviour and creation of a 
corresponding potential well which develops proportion- 
ally to the interaction strength (3 (Fig. 6). At the same 
time the correlation of both phonon modes is involved 
(last term of I©) via the multiple (Rabi) tunneling medi- 
ated by the mode 1 as a source of the nonlinearity. The 
nonlinearity significantly enhances the quantum fluctua- 
tions in the region of maximal tunneling between the levels 
|14| . Competition of the terms responsible for selftrapping 
and phonon-assisted tunneling between the levels result- 
ing in formation of two ground state regions is accompa- 
nied by increased anomalous quantum fluctuations in the 
crossover region. There is to be emphasized that, similarly 
to the case of A = 1 (Fig. 2), there exists no sharp transi- 
tion line in the ground state of the energy in the phase di- 
agram. The transition region is smeared by a width of the 
phonon frequency SI. Sharp transition line occurs rather 
as a well known artefact of some variational approaches 
and of the adiabatic approximation 

In order to describe the quantum fluctuations result- 
ing from the complex interplay of the above described 
contributions we have numerically evaluated the Shannon 
entropies Sq, Sp (Figs. 7a, b) and their sum Sq + Sp 
(Fig. 8a) as functions of fi and x- Here Q = {Qi,Q 2 } 
and P = {-Pi,-P2} and all integrations are meant in the 
two-dimensional space Qi x Q 2 , resp. Pi x P 2 . 

In the weak coupling region, /i = a 2 /2S7 2 < 1, forma- 
tion of the "phases" at % < 1 and x > 1 is greatly reduced 
because of the fluctuations i? and strong correlations be- 
tween the phonon modes (Fig. 8a). For large n the extrema 
at x — 1 g e t sharper but remain non-singular due to the 
finite (although small) quantum fluctuations SI. The obvi- 
ous asymmetry of the sum of entropies Sq + Sp as a func- 
tion of x comes from the fluctuations of the x > 1 " phase" . 
In the limit of large %, both uncertainties, entropic and 
momentum ones, tend to their values corresponding to 
two harmonic oscillators, 2(1 + log7r) = 4.28946 and 0.5, 
respectively (Figs. 8a, b or 9a, c). This resembles much the 
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tunneling dominated region of one-phonon model (Section 
3) with the single exception that the tunneling region is 
spread for all /i because of the pronounced crossover en- 
hanced by phonon-2 pumping. 

Numerical result for the product yj AQ\AP 2 as func- 
tion of (i and x is shown in Fig. 8b and 9b, d. The contour 
plots of Figs. 8 presented on Fig. 10 help visualizing their 
difference, namely (i) the classical contribution of the dis- 
placement growing with the coupling \i to the Heisenberg 
product, Fig. 10. b, and (ii) prevention of formation of 
the sclftrapping dominated and the tunneling dominated 
regions in the weak coupling region \x < 1. Except for 
the region fi > 2 the difference between Fig. 8a and 8b 
is seen not so clear as on the corresponding Fig. 3 a),b) 
for one-phonon model. However it can be easily traced if 
one compares the corresponding crossections along x and 
\i axes, visualised on Fig. 9a)-d). In particular, of main in- 
terest for the application of suggested alternative measure 
are regions of large fj, close to x — 1 an d of x < 1 where the 
wavefunctions exhibit pronounced multipeak structure, as 
seen from Fig. 5a), b). 

The region of maximum in Figs. 8a, b at fi ~ 1 refers 
to the crossover between the quantum fluctuation dom- 
inated region /i < 1 and the selflocalization dominated 
region fi > 1. Namely, at [i ~ 1 there occurs a crossover 
from two correlated oscillators towards two independent 
oscillators, which is visualized in Fig. 8a. These two in- 
dependent oscillators (see wavefunctions on Fig. 5a) refer 
to the complexity of the classical (i.e. adiabatic limit) po- 
tential and the entropic relations handle them accurately 
marking this region as a "classical" one with close to two- 
oscillator value 2(l + log7r). On the contrary, the variances 
in the left-hand side of the Heisenberg uncertainty relation 
(Fig. 4b) continue increasing for small x with growing fj,, 
but apparently weaker than in the one-phonon case of Sec- 
tion 3. In the two-phonon case the separation of two peaks 
of the wave function in the selftrapping region is much less 
pronounced than in the one-phonon case, which is seen 
also from the moderate increase of the product of vari- 
ances at Fig. 8b, 9c, 10b, for x < 1 m comparison with 
the corresponding one-phonon case at Fig. 3b. 

5 Conclusion 

Two-level electron-phonon systems with one (exciton) and 
two (E(g)(6i+&2) Jahn- Teller model) phonon modes exhibit 
reflection symmetry as the source of a hidden nonlinear- 
ity which reveals explicitly by appropriate unitary diago- 
nalization (Fulton-Gouterman transformation ©) of the 
Hamiltonian. This diagonalization (i) excludes electrons 
leaving us with solely phonon Hamiltonian (exact decou- 
pling in electron subspace), (ii) reveals quantum correla- 
tion of phonon modes in the case of two-phonon model. 
As a consequence, related phonon wave functions exhibit 
multipeak, i.e. essentially non-Gaussian structure even in 
the ground state. Besides this "topological" anisotropy, 
additional anisotropy appears in two-phonon model due 
to the correlation of the modes. Resulting complex topol- 
ogy of the ground state implies appearance of two regions 



in the phase plane or /i, [3: the selftrapping dom- 
inated region for x < 1 or 2/i > (3 and the tunneling 
dominated region for x>lor2/i</3 respectively, for 
both models due to competition of two antagonistic (clas- 
sical selftrapping and quantum tunneling inter- 
actions in the Hamiltonian (terms with a and f3 in (0). 
The enhancement of fluctuations close to the crossover 
between these two regions behaves in analogy with the 
corresponding items of the theory of critical phenomena 
but appears smoothed by the finite phonon frequency J?. 
Namely, quantum fluctuations due to finite phonon fre- 
quency Q prevent formation of the ordered "phases" in 
the region of weak coupling jj, < 1 in the phase plane /i, (3 
or /z, x- The interplay of the classical and quantum terms 
yields complex interplay (mixing) of quantum fluctuations 
oc fl and nonlinear fluctuations oc x- 

The moment and entropic uncertainty measures repre- 
sent two classes of characteristics used in the probability 
theory in order to describe quantitatively the spreading 
of the probability distribution for an observable. The un- 
certainty principle in quantum mechanics leans essentially 
on the probabilistic interpretation of the wavefunction and 
can be thus formulated in a twofold fashion - in the form of 
either moment (variance) or entropic uncertainty relation. 

The aim of this paper was to show that for certain 
parameter region the appropriate quantitative measure of 
phonon quantum fluctuations of the systems in consid- 
eration could be the entropic uncertainty measures rather 
than the moment measures (variances) used in the Heisen- 
berg uncertainty relation. We have compared results of nu- 
merical simulations for both types of measures and show 
that the latter exhibits serious shortcomings in description 
of the global fluctuations of such systems. 

Namely, Heisenberg uncertainty relations impose the 
measure under integration (Qf , resp. Pf) which is nonin- 
variant in the space of variables - both translationally and 
rotationally. Hence there are two sources of the disquali- 
fication of the Heisenberg variances, both resulting from 
the anisotropy of wavefunctions caused by the reflection 
symmetry of considered models: (i) former non-invariance 
has its consequences in difficulties describing the multi- 
peak distributions, and (ii) latter one - when trying to 
judge about two or more strongly correlated variables. 

Problem of characterization of quantum fluctuations 
by variances for a system of two correlated oscillators is 
much more subtle than for one oscillator of Section 3. In 
addition to the problems imposed by multipeak distribu- 
tions a serious difficulty arises if several coupled variables 
are involved (in our case - two strongly correlated phonon 
oscillators). If for one- mode case (and for the probabil- 
ity distribution with a single peak) (AX 2 ) could stand 
for an effective "width" of the distribution function, in 
many- variable case it represents only an average "width" 
in the JT-direction, that is the measure under the integral 
is not the best choice since it is strongly affected by, e.g., 
basis rotation. As an example just consider a squeezed 
correlated two-phonon trial function |14| of the general 
form &tp (—^2i j = i 2 a ij x i x j) representing an arbitrarily 
turned ellipsoid in (xi, X2)— plane. To judge about the 
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"width" of such distribution in terms of variances one 
should therefore consider rather tensor quantities repre- 
senting the "width" of the distribution for every direction 

The entropic measure, on the contrary, does not suffer 
from this shortcoming. It presents a handy characteris- 
tics of essentially scalar character which can be applied 
for any number of coupled variables, but is not affected 
by the basis rotation. This holds generally for a measure 
which suggests taking an average of some function of the 
distribution itself (averaging of logP(ir), as for Shannon 
entropy considered here, or generalized entropic measures, 
like, for example, Renyi entropy |21j). 

Similar problems arise, e. g. in quantum optics, at the 
description of fluctuations of the photon coherent Schro- 
dinger cat states and photon multimode correlated sys- 
tems (see, e.g. Ohya et al and Buzek j22| and ref- 
erences therein). The entropic uncertainty relations are 
commonly recognized there as a good tool for handling 
fluctuations of quantum origin. 
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c) u = 2 d) /j = 4 




Fig. 1. Ground state wave functions of one-phonon model parametrized by /3 and /i, /i = 0.6 a); /i = 1 b); (j = 2c);fi = 4 d); 
f? = 1. 
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Fig. 2. Ground state energy in the plane [i and j3. 
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a) S Q + S P c) 12=1 




Fig. 3. Sum of entropies Sq + Sp (a) vs product of variances y / AQ' 2 AP 2 (b) for one-phonon model; Q = 1. Product of 
variances (upper curves) and exponent of the sum of entropies (lower curves) shifted to the respective minimal values 1/2 and 
1 + log7r, respectively, for /3 = 1 (c), f3 = 2 (d). Illustration of the growing product of variances as functions of the classical 
parameter fj, when compared with the sum of entropies. 



yj SP 2 5Q 2 V- 




Fig. 4. Product of variances \J AP 2 AQ 2 calculated from the variational approach of Section 3.1 1121131 . Compared to Fig. 3 
the fluctuations are completely suppressed on the right part of the Figure. 
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a) x=0.9 b) x=l 




c) X=l-5 




4 



Fig. 5. The two-phonon numerical ground state wave functions at fi = 2 and x = 0.9 (a), \ = 1 O 3 ) an d X — 1-5 (c). 




Fig. 6. Effective potential in the plane 71, 72 containing several competing minima for = 2, \ = 1-5- For this values of the 
parameters the narrow (tunneling) minimum 71 ~ 0, 72 > dominates. For \ < 1 the broad minimum at 71 < 0, 72 — would 
dominate (the selftrapping) . 




Fig. 7. Shannon entropies Sq (a) and Sp (b) as functions of x and \i. 
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a) S P +S Q 




Fig. 8. Sum of entropies Sq + Sp in the plane The lower bound corresponds to the two harmonic oscillators value 

2(1 + log7r) (a); Product of variances U AQ'fAPf in the plane x,/- 1 (b). 

The difference of two Figs is not so well pronounced as on its one-dimensional counterpart Fig.3, but for better insight the 
reader is referred to following Fig. 9 with pertaining crossections. 
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Fig. 9. Cross-section of the sum of entropies at Fig. 8a, Sq + Sp as a function of x for various jx (a). The same as a function 
of pi for several x ( D )- Cross-section of the product of variances \J AQ'{APf at Fig. 8b as a function of \ for various n (c). The 
same as a function of fi for various x (d). 
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Fig. 10. Contour plots of the sum of entropies of Fig. 8a (a) and of the product of variances of Fig. 8b (b). 



